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t> ■ 

O ! We study the effects of the finite size of the baryons on the equa- 

tion of state of homogeneous hadronic matter with hyperons. The 
finite extension of hadrons is introduced in order to improve the per- 
fomance of field theoretical models at very high densities, simulating 
the short-range hard-core repulsive part of the baryon-baryon poten- 
tial. We choose the Quark Meson Coupling model to describe the 
baryon dynamics because this model provides a density dependent ra- 
dius of the bag. In addition we calculate finite size corrections for the 
Zimanyi-Moszkowski model in order to make a definite comparison. 
In both cases finite size effects are implemented through a Van der 
Waals like correction in the normalization of the baryon fields. In this 
approach we investigate /3-stable matter in conditions that could be 
found in the interior of massive stars. We find that the excluded vol- 
ume corrections contribute significantly for densities above twice the 
symmetric nuclear matter saturation density, although for the quark 
model the results strongly depends upon the equilibrium conditions 
for the confinement volume immersed in a dense medium. 
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1 Introduction 



Investigation of hadronic matter in extreme conditions of density and tem- 
perature is of relevance in physics since the knowledgwe of its properties will 
eventually lead to clarify some questions about fundamental interactions and 
the recovering of QCD symmetries [[!]. 0- 

The role of strangeness in the evolution of relativistic heavy-ion collisions 
[||] as well as in the dynamics of neutron stars j| has been stated long time 
ago. Near the surface of a neutron star, where the density is of the order 
of the normal nuclear density no, hadronic matter is mainly composed by 
protons and neutrons in (3 equilibrium with electrons, in such a way that 
the whole system remains electrically neutral. At higher densities (2-3 n ) 
hyperons could appear, decreasing the energy per baryon. Further increasing 
of the density could cause the formation of different exotic phases composed 
of mesons condensates, quark matter lumps or quark-gluon plasma. The 
transition among phases are esentially pressure induced, and differs qualita- 
tively when there are one or two conserved charges, e.g. baryon number and 
electric charge ||. 

In most hadronic matter studies, the interacting particles are assumed 
as point-like. In an early work |J, it was claimed the necessity to take care 
of the spatial extension of the hadrons and the influence of this fact on the 
collective phenomena at very high densities. It is worthy to mention that 
there are only a few field theoretical models which consistently includes the 
dynamics of the hadron volume, the most commonly used are the Skyrme 
and bag-like models. In the first case the inclusion of finite baryon density 
effects is not straightforward due to the topological character of its solutions 
[[?]]. On the other hand further refinements of the original MIT bag model 
allows to deal with medium effects upon the hadron structure f5f-fl0|l. Also 



in a recent work |T(J , attempts to include repulsion between overlapping bags 
was done through effective short-range quark-quark correlations. 

Even when finite size hadrons are considered, the associated hard core 
repulsion among baryons in the nuclear medium is frequently not taken into 
account. At densities lower than no finite volume effects are small, but these 
repulsive corrections could be very important at higher densities. For the case 
of hyperon matter, this fact affects the relative population of the different 
baryonic components. 

The main difficulty to deal with a hard core potential is that it does not 
admit a perturbative treatment at high densities. Hence it would be desirable 
to have an effective treatment of these volume effects without leaving the 
simplicity of the quasiparticle approximation. 

We adopt a Van der Waals-like method to take into account finite volume 
corrections. In this scheme the total volume appearing in the thermodynam- 
ical quantities is replaced by the available volume, therefore it has the advan- 
tage that the volume of the particle is the only new parameter introduced to 
represent the hard core repulsive correlations. A similar approach has been 
used in related investigations, as for instance in references [p"Tf-||15|. 

In the present work we want to generalize the excluded volume corrections 
to study the properties of hadron matter with strangeness, and to check 
out to what extent the usual nuclear field effective models can be used to 



extrapolate the hypermatter equation of state at high densities. We introduce 
these corrections at the level of the normalization of the baryon fields. This 
point of view allows us to use the same effective Lagrangians usually applied 
to describe the dynamics of point-like particles. Therefore in the quasi- 
particle scheme of the Mean Field Approximation for homogeneous matter, 
the baryons as a whole will be described by a superposition of relativistic 
plane waves, but they will move only within the accessible volume. 

We shall apply this approach to study the properties of homogeneous hy- 
permatter in /3-equilibrium at zero temperature. Although the formalism is 
stated for the whole baryon octet, in order to clarify the discussions we shall 
only consider A hyperons in our calculations, besides neutrons and protons. 
The equation of state of nuclear matter with lambdas has been previously 
investigated in non-relativistic [O] as well as for relativistic [18| regimes con- 



sidering point-like particles. However a full treatment must include heavier 
hyperons, as deduced for example, from the influence of the H particles on 



the stability of lambdas in medium [19 



In the next section we give a resume of the effective models to be used 
in our calculations. Section 3 is devoted to extend these models to include 
correlations based in excluded volume restrictions. To compare the effects of 
the correction proposed, we have selected two descriptions of very different 
kind. On one hand we select a hadronic Lagrangian for point-like particles, 



the so-called Zimanyi-Moszkowski (ZM) model [2C], for which a fictitious 
size is assigned to the baryons, and on the other hand we have a bag model 
generalized to include quarks interacting with scalar and vector mesons, the 
Quark Meson Coupling (QMC) model |§ §. In the last case two possible 
conditions for the equilibrium of the confinement volume in dense matter are 
proposed. All these items are discussed in section 4. Numerical results and 
discussion are given in section 5, conclusions are drawn in section 6. 



2 Relativistic Effective Lagrangian for Point- 
like Baryons 

Relativistic Hadron Field Theories have proved to be adequate in describing 
many nuclear matter properties at saturation density, and also finite nuclei 
structure, with a relatively small number of free parameters. These theories 
allow to extrapolate matter behaviour to high densities and energies. Al- 
though a chiral treatment is desirable at high densities, we do not consider 
this issue in this work. 

There are a variety of effective relativistic Lagrangians, the more usual 
ones describe nuclear interaction mediated by scalar (a), vector (u) and 



pseudovector (p) mesons |2(J |2T], [22fl . Since neutral matter at high densities 
can be composed not only by nucleons and leptons but also with more massive 
baryons, we select a Lagrangian of the form: 

B 1 A 

+ ~(0"<70„<7 - m a 2 a 2 ) - \F^F^ V + X -m 2 



- lG^G^ + ^m p 2 b^.b M + ^ [¥{iYd,-mi)¥] (2.1) 

In Eq. fl2.lQ F^ v and G M1/ are the tensor strength fields of the uj and p mesons, 
respectively, and the index B(l) runs over all baryons (leptons) considered. 
The term Mb* contains the interaction of the a meson with the baryonic 
fields ty B , and plays the role of the baryon effective mass. The different 
forms of Mb* lead to a family of hadronic Lagrangians. In order to have a 
general discussion Mb* will be kept undetermined until needed. 
We solve the field equations for the case of homogeneous infinite static matter 
in the mean field approximation (MFA) [[H], |22| . In this case all meson fields 
are replaced by their average values, i. e. 
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Up = UoS^o = -^^^g^n B 5^ (2.3) 



w b 

b a , = b 5, 5 a3 = -l^Y.9pfn B (2.4) 
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where a = 1,2,3 runs over all isospin directions and I B is the third isospin 
component of baryon B. In all our calculations we use the values of m a = 
550MeV, = 783MeV, and m p = 770MeV for the meson masses. 
The equation of motion for the baryon field B (spin 1/2) is 

(•Yfy - g B j°u - gflf 7 % - M B *)V B = (2.5) 

The equation above can be used to define the energy /c for a baryon carrying 
momentum k 



k B = ^M B * 2 + k 2 ± g B u ± gflfb (2.6) 

for particle (+) and antiparticle (— ) solutions. Within the MFA at zero 
temperature only the particle solutions contributes. 

The scalar (nf) and baryonic (n B ) densities are defined with respect to the 
ground state of the system \GS > for each baryonic class 



nf =< GS\^ B ^ B \GS >= d -^M B * dk . = = (2.7) 

n Jo ^JMb* 2 + k 2 

k 3 

n B =< GS\& B V B \GS >= -^r (2.8) 
1 1 3tt 2 K ' 

In Eqs. ( |2.7D and ( |2.8| ) ks is the Fermi momentum for baryons and the 
factor *& — 1 for point-like baryons. 

For the Lagrangian (|2.1| ) leptons behaves like free Dirac particles. The 
leptonic density n\ defines the Fermi momentum k\ through 



A-, 3 

n l =< GS\& ¥\GS >= (2.9) 



Mean field equations are solved for fixed total baryon density n and zero 
total electric density charge n e : 



J2n B (2.10) 
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"' ]T<A B + E<? V = (2.11) 
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where g B are the baryon electric charges, q l = — e the electron and muon 
charges. 

All baryon and lepton species are assumed to be in j3 equilibrium. This 
constraint and the condition ( |2.11| ) determine the relative population of the 



particles. In the present work we consider only three baryon species, namely 
B = proton (p), neutron (n), lambda (A), and two lepton species, I = elec- 
tron (e), muon (//). 

Consequently chemical equilibrium is impossed through the following rela- 
tionships among the chemical potentials p: 

= + 

p n = /i A , (if lambdas are present) 

p e = //\ (if muons are present) (2-12) 



with fi B = k§(k B ), see Eq. ([TEfl, and \i l = k l {ki) = \Jmi 2 + k 2 . 

One solves the system of equations ( |2.2|) to ( |2.12|) for a fixed value of the 



total baryon density. This determines the meson fields (ctq, ujq, bo), the baryon 
and lepton densities (n B ,n l ), all the Fermi momenta (ks,ki) and chemical 
potentials (p B , ll 1 ). 

Once the solution has been found, we can calculate the total energy density 

2 2 2 2 2t2 
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and the total pressure 
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1 rffcfc 4 
+ 3^y'o yjmf 2 + fc 2 (2 ' 14) 

3 A Relativistic Effective Approach for Fi- 
nite Size Baryons 



Effective Lagrangians like the one in Eq.( |2.1|) , describe the long range be- 
haviour of nuclear forces. At short distances, i.e. high densities, the strong 



repulsive component of the baryon-baryon interaction appears as a conse- 
quence of the internal structure of the particles. In turn, this short range 
repulsion can be assimilated to a simplified model where baryons are de- 
scribed as extended objects over a spherical volume of radius R, where R 
reflects the range of the core repulsion. Therefore the fraction of available 
space is reduced with respect to the case of point-like particles. A similar 
approach has been applied to study the phase transition of nuclear matter 
to the quark-gluon plasma (JTTJ, [EJ] and in heavy-ion collisions [ ]T3"| , 13 ]. 



In our method we extend the Van der Waals prescription to take into account 
the finite size of different kinds of baryons including hyperons, at any finite 
density. For this purpose we shall use the effective Lagrangian (|2.1|) valid 
for low densities, but the solutions arising from Eq. (|2.5|) will be normalized 
within the available space. Our approach is expected to be valid for densities 
such that the center of mass of baryons be at a distance greater than 2R 
apart. Otherwise the simple quasiparticle picture breaks down. 
Thus we start from the plane wave solutions ip of the Dirac equation (2.5) for 
fermions of spin projection s = ±1/2 carrying momentum k, and normalized 
in a volume V 

^(x) = V-^ 2 u B (k,s)e-^ (3.1) 



where u B (k,s) is the spinor free particle solution in k space [23]. An analo- 
gous expression holds for the antiparticle solution. 

We further assume that if particles have a finite size, then Eq.(3.1) refers 
to the dynamics of its center of mass. Since extended hadrons are assumed 
not to overlap their motion is restricted to the available space V defined as 

|,| 

V' = V-Y,N B v B (3.2) 

B 

with N B the total number of baryons of class B inside the volume V, and 
v b is the effective volume per baryon of class B. Hence if we renormalize the 
particle (antiparticle) wave function replacing V for V, the effective baryon 
fields \I/ read 

V B ( X ) = (V'y 1/2 J2 I a B (k,s)u B (k,s)e~ lk ^ 

k,s 

+ b B \k,s)v B (k,s)e ik ^\ (3.3) 

as a function of the Fock space operators a and b, for particle and antipar- 
ticle respectively. In this way the finite size of the baryons is automatically 
accounted for into the field dynamics. 

It is interesting to note at this point that a more detailed analysis shows 
that when we have a mixture of different baryons, the excluded volume is 
not exactly the same for all species, but depend upon their respective sizes 
[ |Tj| . To simplify the discussions, in this paper we assume that the sizes of 
the baryons are similar for all the species considered (see table 1), and hence 
the excluded volume is, to a good approximation, the same for every particle. 
With respect to the effective volume per baryon v B , this quantity is expected 
to be proportional to the actual baryon volume, i.e. for spherical volumes of 
radius R B 

47T 

v B = a—R B 3 (3.4) 



and for sharp rigid spheres a is a real number ranging from 4, in the low 
density limit, to 3\/2 /ir, which corresponds to the maximum density al- 
lowed for non overlapping spheres, in a face centered cubic arrange. Since 
we wish to apply our method to study the high density behaviour of homo- 
geneous isotropic matter, we shall adopt a = in all our calculations. 
Thus vb = 4v / 2-Rb 3 and the upper limit for the baryon density is given by 
l-AWa; = V2/(8R max 3 ), where R ma x denotes the biggest radius among all 
present baryonic classes. 

In order to see how volume corrections appear in our approach, we shall 
use the renormalized field (3.3) to calculate the relationship between baryon 
densities n B and the Fermi momenta k B 

n B = V~ 1 [ dx 3 <GS\^ B ^ B \GS>=(l-T j n B 'v BI )^4 (3-5) 

JV 37T 

where — > V'/(2tt 3 ) J dk 3 has been used. 

k 

This result is equivalent to Eq. Q2.8Q , if the factor $ takes now the value 

<& = l-Y,n B v B (3.6) 

B 

for the case of finite baryon effective volumes Vb ■ In the limit vb ~^ one 
recovers the point-like expressions. 

Thus we can still use the same expressions we have calculated in section 2.1 
with d satisfying Eq. (|3.6| ) 

Equation ( |3.5| ) shows that volume correlations couples non linearly the baryons 
among themselves, in a density dependent way. 

We can solve explicitly Eq. ( |3.5| ) for n# as a function of all the Fermi mo- 
menta ks, namely 

" B = ~, K < 37) 

(1 + E v B .) 3t 2 

B> 

Since $ in our approach depends explicitly upon the baryonic densities, the 
chemical potentials get an extra term, i.e. 

" B = (£) ^» +A " s (38) 

B'^B 

where fx B = k B (ks), Eq. ( |2.6|) , is the free quasiparticle eigenvalue and 

a b v B ^ [ k B> dkk 4 

3lT b> Jo JM B ,* 2 + k* 



This value for Afi B is in agreement with the one used in the literature |l I 
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The expression for the energy density e is formally the same given in ( 2.13Q , 
and the total pressure p adquires an additional term Ap 

p = p + Ap = p + J2 " B A// B (3-10) 
b 



where po is defined in Eq. Q2.14]) , and is evaluated according to Eq. ( |3.(j| ) 
in all the expressions. 



In our approach only the baryonic states receive an explicit correction due 
to short range forces, through the normalization of the fields \E fB . This can 
be justified because meson fields contribute only through their mean values, 
which are completely determined from the baryonic sources. These sources 
already include finite size corrections in Leptons do not experiment strong 
interactions, and in this sense they are taken as point-like particles. To re- 
sume, in our approach to finite volume correlations, matter properties are 
determined applying the set of equations ( |2.2| ) to Q2.13D , together with ( p.8|) - 
( |3.10|) for a fixed value of the total baryon density n, but using the value of 
■& given in (|3.6|) instead. 



4 Effective Baryonic Models 
4.1 The Zimanyi-Moszkowski Model 

In order to have definite predictions of finite size effects we apply the scheme 
described in section 3, to the ZM model [EDI. This is a local field model 
which basic degrees of freedom are point-like hadrons. This model has been 
extensively used to study nuclear matter and nuclei properties. It has the 
advantage of reproduce the main nuclear matter properties in the lowest 
order of approximation with reasonable accuracy, and using the minimum 
set of free parameters. In the ZM model, properly generalized to include 
several baryonic species, the effective baryon mass Mb* and its derivative 
(1Mb* /da are given by 

M B * = M B (l + g*a/M B y\ 
dM B* R/, , „B 



da 



-g«(l + g »a/M B )- (4.1.1) 



respectively. 

We place this expression for M B * in the equations of sections 2 and 3, to 
calculate the effects of finite size corrections in the thermodinamical bulk 
properties. 

If we take i? = 1, section 2, we describe a system of point-like baryons, and 
using $ from Eq. ( |3.6| ), we include finite size effects. In the last case a ficti- 
tious volume of constant radius is assigned to each kind of baryons in order 
to contrast excluded volume effects with the point-like case. However, the 
high density results must be taken with caution since causality violation is 
expected due to the hard-core repulsion fll3"l . This problem can be avoided 



if the baryon volume decreases with density as predicted by the model pre- 
sented in the next subsection. 
Numerical results will be given in section 5. 



4.2 The Quark Meson Coupling Model 



Another baryon model we shall use is the Quark-Meson Coupling (QMC) 
model H, §]. The QMC model may be viewed as an extension of the Quan- 
tum Hadro-Dynamics (QHD) models like, for instance, the Walecka |?T], E2| 
or the Zimanyi-Moszkowski ^(J models. In the QMC model hadrons are 
represented as non-overlapping bags containing three valence quarks, the bag 
radius changing dinamically with the medium density. Baryons are interact- 
ing by the exchange of a and u mesons that couples directly to the confined 
quarks. It has been found that these extra degrees of freedom provided by 
the internal structure of the baryon lead to quite acceptable values of the nu- 
clear matter compressibility at saturation. Despite the explicit quark fields 
present in the QMC model, hadronic thermodynamical properties are evalu- 
ated in such a way that baryons are handled as point-like particles with an 
effective mass M B which depends on the a field. 

To describe briefly the QMC model, one considers baryons as spherical MIT 
bags where quarks are confined, the Dirac equation for a quark of flavor 
q, (q = u,d,s), of current mass m q and If third isospin component is then 
given by 

(ird„-g^ o co o -g q p Iij%-m q *)W = (4.2.1) 

In this equation all meson fields have been replaced by their respective mean 
field values, in a similar way as we have done in Eq. (|2.5| ) for baryons in 
uniform hadronic matter. Mesons are supposed to couple linearly only to u 
and d quarks, i.e. g s a = g s u = 9 s p — 0. Therefore the effective m q * masses are 

m u / = m Ujd - g u /o 
m s * = m s (4.2.2) 

For a spherically symmetric bag of radius R B , representing a baryon of class 
B, the normalized quark wave functions ty q B (r,t) are given by 

n(r,t) =Af B 1/26 -^- ( - R H X J B ( r/RB) m , ) X* (4.2.3) 



where x 9 is the quark spinor and 



n 



e q = ^ + gl^ + gill b (4.2.4) 



Rb 

Mb = Rb 3 [2n qB (n qB - 1) + R B m q *\ (4.2.5) 

X qB 



0, 



qB 



QqB — Rs^q 



,1/2 



(4.2.6) 



yt qB + RBm q * _ 

with Q q b = [x 2 q B + i.RB m q*) 2 } 1 ^ 2 ■ The eigenvalue x q s is solution of the 
equation 

jo(x qB ) = (3 q ji(x qB ) (4.2.7) 

which arises from the boundary condition at the surface of the bag. 

In this model the ground state bag energy is identified with the baryon mass 

Mb* = + ^ W (4.2.8) 
tin o 



where is the number of quarks of flavor q inside the bag, B^ 4 = 210.85MeV 
is the bag constant and zqb the zero-point motion parameter, fixed to repro- 
duce the baryon spectrum at zero density. 



baryon 


M(MeV/c 2 ) 


R(fm) 


zo 


P 


938.2723 


0.6 


4.00496464 


n 


939.5656 


0.6 


4.00102120 


A 


1115.63 


0.606 


3.86577814 



Table 1: Baryon parameters. M and R are, respectively, the mass and the radius at zero 
density. For the ZM model R is density independent. The bag constant Bo has been fixed 
by the value B 1 ^ = 210.8AiW. 



To obtain this mass one has first to solve Eq. ( 4.2. 7| ), with ( [4.2.2 ) and 



the o"o mean field given by ( |2.2|) where the derivative in this equation must 
be interpreted as {dM B / da) R . 

Equation Q4.2.8| ) shows that the baryon effective mass is a function of the 
bag radius Rb- In the original MIT bag calculations Rb is a constant fixed 
at zero baryon density, but in the QMC it is a variable dynamically adjusted 
to reach the equilibrium of the bag in the dense hadronic medium. We shall 
consider two possibles equilibrium conditions, to be combined with the point- 
like and finite size considerations. They are: 

a) the standard QMC prescription which minimizes Mb* with respect to 
Rb keeping the a field constant, namely [|] 



'dM 



B 



dR f 



(4.2.9) 



This prescription will be denoted as QMCa. 

b) an alternative version of the QMC model, more suitable to account for 
the bag equilibrium at high pressures, is obtained by impossing that the net 
momentum flux through the bag surface be zero (see Appendix) 

-(?££) *L_ (4.2,0) 

where the left hand side is the pressure exerted by the quarks inside the bag 
of kind B and the right hand side is the baryon contribution to the total 
hadronic pressure. 

This condition will be labelled as QMCh in what follows. 

It must be mentioned that a similar condition can be obtained by min- 
imizing the energy density e with respect to the effective bag volume vb, 
keeping constant the total volume V, namely |l~2]j 



nffdM B *\ 1 ^ r h B> dkk 



— V / . (4.2.11) 
3tH - Jo Jm* 2 j_ w 



The only difference with Eq.( |4.2.IUD stands in the factor nf /n B . This factor 
is found numerically to be close to unity in the range of densities we handle 



in the present work. In fact, numerical results using condition ( [4.2.11 ) dif- 
fer only by negligible amounts with respect to the predictions arising from 
([4.2.10|) for symmetric nuclear matter. 

Hence we prefer to use the much simpler condition ( |4.2.10| ) for further cal- 
culations. 

Conditions QMCa and QMCb coincide only in the case of vanishing 
density, but the last one seems to be more appropriate to describe general 
equilibrium situations of bags inmersed in high density medium. 



Once Mb* has been defined microscopically, the dynamics of the hadrons 
in the QMC model arises from the effective Lagrangian ( |2.1| ), as discussed in 
section 2. 

Only the couplings of mesons to baryons remains to be properly defined, since 
in this model mesons interacts directly with quarks. In fact, the quark-meson 
couplings g^, with = a,u, p and q — u,d, are related to the corresponding 
meson-baryon couplings in a simple way || . If we denote with g B the coupling 
of the 0-meson to 5-baryon, and calling n B s the non-strange quark content 
of baryon B, then 



9(T 


ns u 

3 y ° 


9* 


""ns u 
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= 9 U P 



When $ = 1 (see section 2), we neglect the particle size in the evaluation 
of many-body effects. Otherwise these are included using Eq. ( |3.6| ) for $ 
instead. 

An interesting fact related to the QMC model is that the bag volume reflects 
the changes in the medium density through the mean value of the a field 
appearing in either of the equilibrium conditions QMCa or QMCb. This is 
a self-consistency condition since a itself depends on the effective volume of 
the baryons through the coefficient i?. 
Numerical results will be discussed in the next section. 



Model 


Case 


9a 


9u 


9 P 


9a 


9% 


ZM 


NC 


7.8449 


6.6710 


9.4802 


4.6760 


4.9620 


vc 


7.4039 


5.8275 


9.2446 


4.1997 


4.2872 


QMCa 


NC 


17.964 


9.012 


9.220 


11.977 


6.008 


VC 


17.007 


8.172 


8.873 


11.338 


5.448 


QMCb 


NC 


17.985 


9.024 


9.219 


11.990 


6.016 


VC 


17.046 


8.199 


8.870 


11.364 


5.466 



Table 2: Baryon-meson couplings ga,u>,p used in our calculations, for the ZM, QMCa and 
QMCb models as described in the text, and combined with the point-like (NC) and finite 
size (VC) cases. The couplings has been adjusted in order to reproduce the symmetric 
nuclear matter properties at saturation in either case, see section 5.1. 



5 Finite size effects: Numerical Results 



5.1 Homogeneous symmetric nuclear matter 



Nuclear matter properties has been extensively studied using point-like hadronic 
models. We wish to investigate excluded volume corrections in the homoge- 
neous symmetric nuclear matter equation of state (EOS). With this purpose 
we shall use alternatively the hadronic models ZM and QMC to apply the 
scheme described in section 3. Physical predictions for point-like (i? = 1) 
and finite extension baryons will be compared, and for the QMC case we 
also have considered two alternative conditions for the equilibrium of bags 
in the nuclear medium. In the following we shall use the labels VC and 
NC to indicate results obtained with or without excluded volume correction, 
respectively. Also the labels a and b will be associated with the QMC model 
to distinguish between the proposed equilibrium conditions, as explained in 
the previous section. 




Figure 1: Thermodynamical properties of symmetric nuclear matter: binding energy 
Eb = (e/n) — Mc 2 (a) and pressure (b) as functions of the relative baryon density. The 
curves corresponds to the models Zimanyi-Moszkowski (ZM) and Quark-Meson Coupling 
(QMC). Results with excluded volume correction (VC) are represented with dashed lines 
and those without this correction (NC) by solid lines. For the QMC model the cases 
QMC a and QMCb arise by selecting different equilibrium conditions as explained in the 
text. 



The equations to be solved for nuclear matter are the same as those in the 
neutral matter case described in section 2, except that proton and neutron 
densities are equal (n p = n n ), hyperons and leptons are absent and hence 
Eq. (|2.11 ) for electrical neutrality no longer holds. 



In tables 1 and 2 we display the set of parameters to be used in both the 
ZM and QMC models, respectively, combined with different possibilities, as 
explained in the table captions. Each set of parameters has been obtained 
with the constraint that the known symmetric nuclear matter properties at 
saturation must be reproduced, i.e. the saturation density no, the binding 
energy per nucleon E b0 , and the symmetry energy coefficient a s : 

n = 0.15/m~ 3 

E b0 = (e/n) - Mc 2 = -WMeV 



- 1 ^L- w " v (5 ' Ll) 

where n = (n n + n p ) is the total nucleoli density, t = (n n — n p )/n and 
M = 938.92MeV/c 2 is the average free nucleon rest mass. 
When finite size baryons are considered we fix the nucleon radius in vaccum 
at 0.6/m; with this choise the upper nuclear density allowed for non over- 
lapping bags is always beyond n/n = 5, for the set of parameters in tables 
1 and 2. 

In figure la we plot the binding energy Eb = (e/n) — Mc 2 , as function of the 
relative baryon density n/n up to 5 times the saturation density, for both 
ZM and QMC models, with and without finite volume corrections. Also we 
show in figure lb the total pressure for symmetric nuclear matter, corre- 
sponding to the same cases depicted in figure la. We can see in either figures 
la and lb that all the cases practically coincide in the range < u/uq < 1.5. 
Thus, as expected, in this density domain the short range nucleon-nucleon 
repulsion can be neglected. 

In particular the nuclear matter compressibility at saturation, k = 9(dp/ dn) no , 
is only slightly increased by finite volume corrections, the increment is about 
15% in the ZM and about 8% in the QMC models, see table 3. 



Model 


Case 


M p */M p 


k (MeV) 


ZM 


NC 


0.850 


221.4 


vc 


0.862 


254.1 


QMCa 


NC 


0.773 


300.6 


VC 


0.792 


323.4 


QMCb 


NC 


0.772 


296.1 


VC 


0.791 


318.6 



Table 3: Predicted properties for symmetric nuclear matter at saturation, using the ZM 
and QMC models with the corresponding set of couplings constants detailed in Table 2. 

For the QMC formalism the cases a and b are practically undistinguishable 
when volume corrections are not included, meanwhile in the VC approach 
there are significative differences only for densities higher than 2.5n . For 
every case the inclusion of the proposed corrections enhances, for a given 
density, both pressure and binding energy. This increment is greater for 
QMCa than for QMCb. This behaviour can be explained by looking at fig- 
ure 2a, where we plot the density dependence of the proton radius R p . In 
fact, one sees there that all predictions are very similar for densities below 
no, but at higher densities the bag size decreases much faster in the option 
QMCb including or not volume corrections, as compared with the option 
QMCa and therefore in the option QMCb d is always closer to unity. The 
behaviour shown in figure 2a could be expected because in the condition 
QMCb the bag must support the growing external baryon pressure, while 
option QMCa only equilibrates the in-medium bag with its internal struc- 
ture. 

To conclude, in case QMCb the excluded volume remains negligible in a 
wider density range, allowing the baryonic system to compress to much higher 



densities as compared with the standard option QMCa, or with the rigid, 
undeformable baryon used for the ZM model. 

Also it must be noted in this figure that corrections effects have opposite 
signs in cases QMCa and QMCb. 

The pronounced compression of the bags predicted by the QMCb case could 
justify the use of effective point-like equations to extrapolate thermodynam- 
ical bulk properties at high densities, provided the bag compression is appro- 
priately accounted for in the effective baryon mass. 
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Figure 2: Properties of nucleon immersed in symmetric nuclear matter: the proton 
radius (a) and the proton mass (b) as functions of the relative baryon density. In the 
first case only QMC results are shown since in the modified ZM model a constant radius 
is assumed for all the baryons. The lines and labels conventions are the same as those 
defined in fig. 1. 

In figure 2b we plot the behaviour of the proton effective mass M*/M p as 
function of the relative nucleon density. This figure shows that repulsive finite 
size effects tends to stabilize the effective mass when the density increases, 
for both the ZM and QMC models. 

On the other hand the difference between equilibrium conditions QMCa and 
QMCb becomes significative for the effective mass only for extremely high 
densities. 



5.2 Homogeneous neutral matter 

We proceed to calculate numerically the effects of finite baryonic size on neu- 
tral matter properties, using alternatively the ZM and QMC models. The 
parameters inherent to each model have been fixed in the previous section, 
tables 1 and 2. 



To simplify the discussion, we study hyperon neutral matter compossed 
only by protons, neutrons and lambdas, together with electrons and muons. 
In the ZM case the coupling of the hyperon A to mesons must be determined 
from experimental data. We use the condition that the binding energy of a 



single A in symmetric nuclear matter must have a minimum of —28MeV |25 



at the saturation density n Q . This condition determines the values g£ , which 



1S00 




Figure 3: Thcrmodynamical properties of electrically neutral nuclear matter with A 
hyperons in (3 equilibrium with leptons: the energy per particle (a) and the pressure (b) 
as functions of the relative baryon density. The lines and labels conventions are the same 
as those in fig. 1. 

are quoted in table 2, for the NC and VC cases. In the QMC framework 
the A coupling constants are reduced to 2/3 of the value corresponding to 
non-strange nucleons; numerical values of g£ w are quoted in table 2. 

In figures 3a and 3b the relativistic energy per baryon e/n and the pres- 
sure p are drawn versus the relative baryon density, respectively. Most of the 
general conclusions given for nuclear matter can be repeated here: the VC 
approach increases, for a given density both pressure and energy, indepen- 
dently of the model used. 
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n/n n/n 

Figure 4: Properties of the A in neutral nuclear matter in (3 equilibrium: the bag radius 
(a) and the effective mass (b) are shown as functions of the relative baryon density. In the 
first case only QMC results are displayed. The lines and labels conventions are the same 
as those in fig. 1. 

The effects of this correction are more significative for ZM and QMC a 
and less important for QMCb. 

The density dependence of the raddii and effective masses of nucleons are 



very similar to the result obtained for symmetric nuclear matter. The be- 
haviour of Ra with the density is depicted in figure 4a for the QMC model, 
and the density dependence of the A mass appears in figure 4b; it can be 
seen that finite volume corrections enhance this effective mass for all the 
cases considered, stabilizing its high density behaviour. 
An interesting aspect of the volume effects arises in the relative baryonic pop- 
ulations, n B /n, and leptonic populations, rii/n, as a function of the baryon 
density, figures 5a and 5b. In figure 5a we study the change in the relative 
particle populations for the ZM model, one can see that the finite volume 
correlations enhance the production of hyperons as compared with the point- 
like case. Furthermore the A particles appear at a relative density of 2.5 (thin 
medium-dashed line), a bit earlier than in the point-like description (thick 
medium-dashed line). This could be expected since at high densities the 
finite size correlations favours the creation of heavier electrically neutral par- 
ticles, decreasing the pressure of the system. This last statement proves to 
be true provided the size of the hyperon is similar to the size of the nucleons, 
as occurs in the present case. 

In figure 5b we plot the relative particle population for the QMCb treatment, 
a qualitative agreement with this figure is obtained when the stability con- 
dition QMC a is used. It can be seen that in the QMC model the inclusion 
of finite volume correlations modifies very little the relative presence of the 
hyperon A, which appears at lower densities with respect to the predictions 
of the ZM model. 

These results can be also compared with those obtained in ref.[[2(J, where 
baryons are considered as point-like. 
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Figure 5: Relative populations of baryons and leptons as functions of the relative baryon 
density for the modified ZM model (a) and the QMC model (b). In the last case only 
results for the equilibrium condition QMCb are shown. Different line types are used for 
each particle. For each class of particle thick (thin) lines represents NC (VC) calculations. 



6 Conclusions 



We have introduced finite size volume correlations to extend the study of 
matter properties at densities much higher than that of symmetric nuclear 



matter at the saturation point. This has been performed through a Van der 
Waals approach to describe hypermatter at zero temperature. One of the 
advantages of the present approach consists in that we can continue to use 
the relativistic effective Lagrangians appropriate for low baryonic density, 
because the excluded volume corrections appear within the normalization of 
the baryon fields. 

We have choosen two effective models to compare the effects of finite volume 
correlations with respect to the more usual point-like descriptions, these are 
the ZM and the QMC models. First we have enlarged the range of applica- 
bility of Van der Waals volume correction to higher densities by introducing 



a model (QMC) that avoids causality violation [13]- [15 1, since the hard core 



has been replaced with a compressible bag structure. Second, although the 
QMC model deals with finite volume bags, this property is usually lost in the 
calculation of thermodynamical properties. This feature has been partially 
corrected in the results QMC-VC of this work. 

Near the normal saturation density finite volume effects prove to be small 
for both the ZM and the QMC models. At higher densities these repulsive 
short range effects show up in a model dependent way. For the ZM model 
with rigid hadrons, the pressure and the energy raise steeper than for the 
point-like prediction as the density gets closer to the touching limit. This 
favours the production of A hyperons in our example. 

For the QMC model we have obtained that the contribution of finite effects 
depends significantly on the bag equilibrium condition used, as explained in 
the text. Using the standard condition the bags are not seriously compressed 
with density and volume corrections can be appreciable in the EOS but not 
in the change of relative particle occupation. 

On the other side, if the bags must equilibrate the external baryonic pres- 
sure, a seemingly more physical condition, their respective radius suffer a 
faster reduction when the density increases. Therefore assuming the validity 
of this approach, short range correlations could be neglected in the EOS up 
to much higher densities, provided the compression effects are included in 
the baryonic effective mass as the QMC model does. In fact those physical 
aspects related to the size of the hadrons itself will be seriously affected by 
the choise of the equilibrium conditions. 

A general conclusion is that finite volume corrections are less significative 
in the EOS when the medium density influence on the hadron structure is 
properly taken into account inside the model prescriptions. 

The production of hyperons is sensibly affected by excluded volume corre- 
lations, decreasing its density threshold and enhancing its relative population 
as density increases. 

This studies on finite baryonic size effects on hyperon matter properties 
aim to describe more properly the possible composition and structure of 
neutron stars. However a more complete treatment must include the full 
baryon octet. 



Appendix 



In order to deduce the equilibrium condition proposed for case QMCb we 
consider only a and uo mesons. The inclusion of other mesons is straigth- 
forward. The QMC lagrangian density with quark fields ^ q {x) coupled to 
scalar a(x) and vector cu^x) neutral mesons, is written as follows 

C-qmc(x) = Ej=i £q(x) + C° Mesons (x) , (A.l) 
C q (x) = (C° q (x) -B)G v - ±%(x)%(x)A s , (A.2) 

Here Gy is the radial step function which schematically confine the quarks 
inside the spherical bag. The terms proportional to the surface delta function 
As ensures a zero net flux of quark current through the bag surface. In this 
lagrangian we have defined the following terms: 

C° q (x) = <I> q (x)(ird,-m q + gla(x)-g«r^(x))* q (x), 
C° Mm (x) = \[d^{x)d,a{x)-mla\x)\-\F^{x)F lxv {x) 

From this lagrangian the equations of motion of quarks and mesons can be 
derived 

9 v (i @- mq + g> - g« </j)y q {x) = 0, (A.3) 

(□ + m> = Ji lS «M^ ( A -4) 
dpFv + m 2 u u" = EU g«%r% ®v, (A.5) 

these equations must be completed with the boundary conditions at the bag 
surface for the quark fields : 27 • n\I/ ? (r = R) = ^ q (r = R) and i^f q (r = R) = 
$> q (r = R)j ■ n. 

In the following equations summation over quark-flavor indices is as- 
sumed. 

The energy- momentum tensor is evaluated by the canonical procedure 
giving 

t qmc = \iQvi^-f d v *) + d»od v o - F<* x d u Lj x - g^C (A.6) 
inserting the explicit form of C and by taking its divergence one obtains 

diFfrc = jiQvd^rf F %) + \in,A s ^ q r $ %) + ±F(* q %A 8 ) 
-e v d^ q tf, + ^ q (-m q + glo - gl 0)* q ] 
+Bn»A s - [^W q ^ q + + &o - gl ^ q ]n v A s 

-d^ x d u u x + Dad" a + \d u [m 2 y + m> A <9^ A ] , (A.7) 



Using the equations of motion and rearranging terms one gets 



^Tqmc = -6v* q d v (g«(T - g* <f>)V + \iA s n^ ql » d v *) + ±d»(%%A s ) 

-[\M q ?V q + V q (-m q + g%a - g* 0)V q ]ri , A 8 
+Bn u A s + (□ + ml)a d v o - (dfiF^ + mluj x )d v u x (A. 

the two terms in the squared bracket are canceled by using (|A.3|) . If we 
denote by Qy the radial step function complementary to Qy, introducing 
the separation 1 = By + Qy in the last line and reagrouping terms 

-Q v [(d/iF^ + mlu x ) - ^g^gjd^ + Bn"A s 

+^A s [in^ q Y d v * fl )] + q %A s ) 

+ [(□ + ml)a d u a - {dfiF^ x + mlu x )d v u x )Q v . (A.9) 

The first two terms on the right hand side are identically zero by virtue of 
Eqs. (|A.4 ) and ( |A.5[ ) respectively, therefore it becomes 



d,T^ MC = A s (Brf + Un^ q rd^ q ) 

+±d"(* q * q A s ) + d,T^ esms Q v (A.10) 

where 

d,UL ons = (D + m^^-^ + my)^ (A.ll) 
In this way we can rewrite the above expression as follows 

n ^t^QMC = n vdpiTbag + n vd^ Mesons® V (A. 12) 

with 

nudX^ = {-P D + B)A s + l -n-d^ q m q A s ) (A.13) 
where P D = -l/2n • d(^ q ^ q )\ 

surface is the pressure exerted by the quarks on 
the bag surface. The standard QMC equilibrium condition can be derived 
by requiring to be zero the first term of Eq.( [A.12| ). Particularly, for the 
QMC model the second term of this equation is strictly zero, which is due 
to the fact that in the QMC lagrangian there is a static bag surrounded by 
free mesons. As the baryon density increases this picture is not adequate, 
since at high densities the mesons are far from the free field approximation. 
Although mesons are in the whole space, i.e. outside as well as inside the 



bags, the Qy in Eq. (|A.12 ) is coming from a cancellation of identical terms as 
it has been shown. The fact that mesons are also inside the bag is reflected 
on the baryon density dependence of m q , x q , R, the quark wave functions, 
etc. 



On the other hand, from the effective Lagrangian [TT] we have Tf v 



hadron 



Tj^ + Tq U , where comes from the free meson sector and T$ v is the con- 
tribution from baryons coupled to mesons. Energy-momentum conservation 



implies d^T^ dron = 0. Explicit evaluation shows that d^T^ = d^T^ esons . 
Consequently, in this limit the second term of Eq.( |A.l^ ) can be replaced by 



—d^T^. Multiplying the resulting equation by n v and integrating on the 
external bag volume we obtain 

4ttR 2 (~P d + B)= J dSn^n u Tg u . (A.14) 

For homogeneous hadronic matter Tq V can be assumed as constant on the 
bag surface, obtaining the following equation 

- P D + B = n^T^r = R) (A.15) 

where n\ is a (space-like) radial unit vector. Therefore, we get 

-Pd + B = W T o J (A.16) 

In the case of homogeneous static matter the space sector of the relativistic 



energy-momentum tensor is diagonal ||27|| , which implies that 



1 3 



-P D + B = -J2T™ = -P (A.17) 

6 i 

where Po is a pressure related to T . In this form we have obtained the 
equilibrium condition of case QMCb. 
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